Inelastic and elastic collision rates for triplet states of ultracold strontium 
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We report measurement of the inelastic and elastic collision rates for *®Sr atoms in the {5s5p)^Po 
state in a crossed-beam optical dipole trap. This is the first measurement of ultracold collision 
properties of a ^Pq level in an alkaline-earth atom or atom with similar electronic structure. Since 
the (5s5p)^Po state is the lowest level of the triplet manifold, large loss rates indicate the importance 
of principle-quantum-number-changing collisions at short range. We also provide an estimate of the 
collisional loss rates for the (5s5p)^P2 state. Large loss rate coefficients for both states indicate that 
evaporative cooling towards quantum degeneracy in these systems is unlikely to be successful. 

PACS numbers: 34.50.Cx 



Metstable ^Pj states of alkaline-earth atoms and atoms 
with similar electronic structure (Fig. [1]) display vastly 
different optical and ultracold collisional properties com- 
pared to states found in alkali-metal atoms more com- 
monly used in ultracold atomic physics experiments. The 
extremely long-lived ^Pq states in Sr and Yb serve as 
the upper levels in neutral-atom optical frequency stan- 
dards [11]. The weakly allowed ^Sq-^Pi intercombination 
transition serves as the basis for powerful laser-cooling 
techniques '3] and may enable useful optical tuning of 
the ground-state scattering length '^P2 atoms inter- 
act through long-range anisotropic interactions that 
allow magnetic tuning of the interactions and cause rapid 
inelastic collisional losses [1, HQ- states of alkaline- 
earth atoms have also been proposed for lattice-based 
quantum computing . 

Here, we report measurement of the inelastic and elas- 
tic collision rates for *^Sr atoms in the {5s5p)^Po state 
in a crossed-beam optical dipole trap. This is the first 
measurement of the ultracold collisional properties of a 
■^Po state, which is of great interest because of its role in 
optical clocks. We also report an estimate of the colli- 
sional loss rates for the *^Sr (5s5p)^P2 state. Large loss- 
rate coefficients for both states make efficient evaporative 
cooling in the ultracold regime unlikely. 

Early laser-cooling experiments with Sr [T]| and Ca 
[T^ l showed that it is straightforward to magnetically 
trap metastable ^P2 atoms through natural decay in a 
magneto-optical trap (MOT). This generated interest in 
the possibility of achieving quantum degeneracy in this 
state and motivated calculations that found novel colli- 
sional properties of the metastable ^Pj levels Mag- 
netic dipole-dipole and electric quadrupole-quadrupole 
interactions between ^P2 atoms produce anisotropic, 
long-range potentials with bound states and collisional 
rates that can be tuned with magnetic field. Reference 
showed that s-wave colliding states can be coupled to 
much higher partial waves of outgoing channels of other 
magnetic sublevels even if the initial state is maximally 
spin polarized. This leads to two-body inelastic loss rates 
in magnetically trapped samples ,6] that are on the or- 
der of elastic collision rates, making efficient evaporative 
cooling towards quantum degeneracy of ^P2 atoms in a 
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FIG. 1: Strontium atomic levels. Decay rates (s~^) and exci- 
tation wavelengths are given for selected transitions. Laser 
light used for the experiment is indicated by solid lines. 
Dashed lines indicate transitions that are only made due to 
spontaneous decay. All levels with energy below 43, 000 K are 
shown. The energy of the ^Si state is 41780.5 isT. The inset 
gives the energy splittings of the metastable ^Pj levels. 



magnetic trap unlikely. This was confirmed in experi- 
ments with magnetically trapped Ca . 

Theory [1 i, only considered Zeeman-sublevel- 
changing (ZSLC) collisions mediated by long range in- 
teractions, which led to speculation ^ 0| that losses 
in optical traps might be low enough to allow evapora- 
tion. However, large loss rates were also found in gases 
of ^P2 Yb atoms held in an optical dipole trap, which 
suggested that fine-structure-chang^ing (ESC) collisions 
at short-range are also significant ^8|. 

By studying ultracold collisions between ^Pq atoms, 
which occupy the lowest level of the fine-structure triplet, 
we remove the possibility of ZSLC and ESC collisions 
and probe the importance of principle-quantum-number- 
changing (PQNC) collisions. Eor ^Pq atoms, PQNC colli- 
sions (also known as energy-pooling collisions [l3|) result 
in one atom in the ground state. Differences between 
■^Po and ■^P2 collisions may also arise because the ^Pq 
atom is isotropic and lacks magnetic-dipole and electric- 
quadrupole moments. ^^Sr has nuclear spin 1 = 0. 
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FIG. 2: Timing diagram for measuring metastable collision 
dynamics (color online). Note the breaks in the time axis. Si- 
multaneous application of the 688 nm and 689 nm lasers popu- 
lates the metastable triplet levels. A strong 461 nm beam ap- 
plied after the initial optical pumping ensures that no atoms 
remain in the ^So state. For measuring ^Po dynamics and to 
keep the level unpopulated, the 3 /xm laser is on the entire 
time, and the 679 nm is pulsed on to repump the ^Po atoms 
to the ground state. For measuring dynamics, the timing 
of these two lasers is interchanged as indicated in the timing 
labeling. 



At higher temperatures, de-excitation of Sr(5s5p)^P/ 
states due to collisions with background noble gas atoms 
and ground state Sr has been well-studied. 
Energy-pooling {5s5p)^Pj + {5s5p)^Pj> (55^)^5' + 
{5s6s)^'^S collision rates have also been measured in a 
Sr heat pipe All of these studies, as well as ultra- 
cold experiments with levels in other atoms 0, S], 
worked in a regime where fine-structure changing colli- 
sions were energetically allowed. To our knowledge, we 
report on the first study that isolates ■^Pq — ^^o colhsions. 

Laser-cooling and trap ping aspects of the experiment 
are described in fll|, Il7j . To obtain higher density and 
longer sample lifetime, and to trap all electronic states, 
^5*0 atoms are loaded into an optical dipole trap (ODT) 
of controllable trapdepth as described in jT3|, yielding 
equilibrium temperatures of between of 3 and 15 /xK and 
peak densities up to 10^'* cm"'^. The potential seen by 
the atoms is determined from measured ODT laser prop- 
erties, well-known atom polarizabilities [U, U^, and mea- 
sured trap oscillation frequencies [l9l |. 

We then pump atoms via ^Sq ^Pi '^Si followed 
by natural decay to the metastable state of interest by 
applying a 689 nm laser for 3 ms while a 688 nm laser is 
on (Figs.[T]and[2|). Atoms decaying from the ^Si state to 
the wrong metastable state are repumped with a clean- 
up laser at 3.0 /im for experiments with ■^Pq atoms and 
with a 679 nm laser for experiments with ■^P2 atoms. The 
clean-up repumper is turned on during the loading stage 
of the ODT and left on for the rest of the experiment. 
Any atoms remaining in the ground state are removed 
from the trap with a 2 ms 461 nm pulse. We typically ob- 
tain 10^ metastable atoms at a temperature near 10 /xK, 



a density as high as lO'^'^cm"'^, and phase-space density 
as high as lO^'^. 

The zero of time for interaction studies is set after the 
metastable atoms have equilibrated for approximately 10 
ms after the 689 nm laser is extinguished. The end of 
the interaction time is determined by the extinction of 
the ODT, after which the atoms ballistically expand and 
fall under the influence of gravity. The density drops 
rapidly enough that atom-atom interactions cease on a 
millisecond time scale. 

The number of atoms and sample temperature are 
determined with time-of-flight absorption imaging using 
the ^Sq-^Pi transition. This necessitates repumping the 
atoms to the ground state during the ballistic expansion 
by applying the 679 nm laser in the case of ^Pq atoms 
and 3.0 laser in the case of ^P2 atoms. It is im- 
portant to release the ODT before repumping because 
the repump lasers cause density-dependent light-assisted 
losses. Complete repumping requires a few milliseconds. 

Photon recoil during optical pumping to the 
metastable state before the interaction time and then 
to the ground state for imaging affects the atomic mo- 
mentum distribution, which complicates measurement of 
the temperature. The recoil energy for a red photon is 
h^k"^ / {mkB) = 0.5 /zK, and from the branching ratios 
shown in Fig. [1] each atom is expected to scatter about 
5 red photons during initial excitation and 1-3 photons 
during repumping. The repump lasers are all aligned hor- 
izontally, so we expect several /iK of extra heating along 
this axis. Recoil from the 3 /.tm photons is insigniflcant. 

An approximate description of the loss of atoms during 
the interaction time can be derived from a local equation 
for the evolution of the atomic density, fi — —(in? — Fn, 
which, assuming constant sample temperature, can be in- 
tegrated spatially to yield the evolution in atom number 
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A'o is the number at the beginning of the interaction time, 
and the one-body loss rate, F is due to background colli- 
sions. The effective volumes weighting one- and two-body 
processes, indexed by subscript q, are defined by 



V,{T)= J dV[n(r)/no]^ 



(2) 



where n(r) is the spatial density distribution calculated 
from the trap potential and atom temperature including 
effects due to truncation of the Boltzmann distribution 
[2l| . and no is the peak density in the trap. The spatial 
integral extends over the region contained in the trap. 

Figure [3] shows representative data for temperature 
and number evolution for atoms in the '^Pq state. The 
horizontal axis {x) is about 2 /iK hotter than the verti- 
cal, which we attribute to photon recoil heating during 
optical pumping. The temperature change during the in- 
teraction time and the x — y temperature difference are 
both small, so an assumption of constant temperature 
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FIG. 3: Left: Number of trapped atoms as a function of 
time for studies of the "^Po state (color online). Two-body 
fits allowing for one and two-body losses using Eq. [T] and 
the numerical simulation described in the text are identi- 
cal. The statistical uncertainties in the two-body loss rate 
constant, /? = 18 ± 3 x 10~^* m'^/s, and one-body loss rate, 
r = 1.4 ± 0.9 s~^ are dominated by correlation between the 
parameters, but /3 is determined much better than systematic 
uncertainties described in the text and is the dominant loss 
mechanism. The one-body fit sets 13 = Q and is not a good 
description of the data. The observed value for T is consis- 
tent with the decay rate for ^5*0 atoms, which is dominated by 
background gas collisions. Right: Temperature evolution of 
the trapped sample. The fit is from the numerical simulation. 
The ODT trap depth resulting from symmetric saddle-points 
located approximately along the horizontal axes (a; and z) is 
Utrap/ks = 20 /iK. The potential barrier for escape along the 
vertical (y) direction is SOfjK. 
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FIG. 4: Two-body collision rate constant measurements for 
the ^Po state (color online). The solid line is P^i = o-iJ\/2, 
where a — Sna^ , \a\ = 100 ao, and v = •^/S/csT/Trm. 



is a reasonable approximation, and Eq. [T] fits the data 
well. The sample temperature used to calculate effec- 
tive volumes is taken as the average of Tx and Ty. A fit 
neglecting two-body decay reproduces the data poorly, 
confirming the importance of two-body effects. 

Values of /3 derived from a series of decay curves are 
shown in Fig. HI The dominant uncertainties in /3 are 
systematic and reflected in the error bars. Uncertainty in 
temperature of 2 //K only contributes a 10% uncertainty 
in /3, while uncertainty in the trap oscillation frequencies 
at the 20% level causes a 60% uncertainty. At our level of 
uncertainty, we observe no temperature dependence and 
find P = (2.3 ± 1.4) X 10-1^ m^/s. 



As pointed out in [20|, both inelastic-coUisional loss 
and evaporation due to elastic collisions are two-body 
loss processes contributing to f3 — (iin + fPei, where 
f3in is the inelastic-collision-loss rate constant, f3ei is the 
elastic-collision rate constant, and / is the fraction of 
elastic collisions resulting in an evaporated atom. The 
rate constants can be expressed as = /3/(l -I- /7) and 
(3ei = 7/3/(1 -|- /7), where the ratio of rate constants is 
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E is the average energy of the trapped atoms, Ei„ is 
the average energy of atoms lost due to inelastic colli- 
sions, and Ef^i is the average energy of atoms lost due 
to evaporation. Parameters on the r.h.s of Eq. [3] can 
be calculated from the trapping potential and equilib- 
rium atom temperature. In our case, these quantities 
must be evaluated numerically because the potential is 
not amenable to analytic solutions and truncation of the 
Boltzmann distribution is significant [21]. For the var- 
ious conditions studied here the ranges of parameters 
are 0.7 < / < 0.3 and 0.6 < 7 < 3.7, and we obtain 
Pin = (1.9 ± 1.2) X 10-1^ m3/s (Fig. HI). The data sug- 
gests that Pel increases slightly with temperature, as one 
would expect since P^i = a-v\/2, where a is the elastic 
scattering cross section and v is the mean thermal veloc- 
ity. The error bars in these measurements are dominated 
by the trap uncertainties. Our treatment assumes er- 
godicty, which is reasonable given the small difference 
between and Ty. 

To check that the assumption of an equilibrium tem- 
perature does not introduce significant bias, we adapted a 
numerical model of the dynamics of atoms in an ODT (2^ 
designed for a system with high 77 = Utrap/ksT to our 
low-ry situation by incorporating the evaporation treat- 
ment of ^211 and numerically calculating all effective vol- 
umes and related integrals using the actual trap shape 
and a truncated momentum distribution. Fits to both 
number and temperature evolution (Fig. [3]) yielded the 
same values of Pin and Pel as above. Assuming pure s- 
wave elastic collisions, this simulation also allows us to 
provide an estimate of the magnitude of the •^Pq scatter- 
ing length \a\ = 100 ± 50 Oq. 

Measurement of the number and temperature evolu- 
tion for atoms in the state, together with numer- 
ical simulations assuming an equal trap depth for all 
magnetic sublevels, constrain Pin — (1.3 ± 0.7 ± 0.7) x 
10~^^ m^/s reasonably well. The first uncertainty comes 
from the fit, and the second uncertainty is systematic 
from knowledge of the trap potential. Several factors 
that increase uncertainty in determination of P in this 
case must be addressed. The temperature measurement 
along X is distorted by a high density of atoms evaporat- 
ing into the arms of the crossed-dipole trap. However, 
the numerical simulation and elastic-collision rate im- 
plied by the low Ty indicate that the temperature change 
and lack of equilibrium is not significantly worse than 
in ^Pq. The different magnetic sublevels have different 
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AC stark shifts and experience optical trap depths that 
vary by almost a factor of two. This complicates model- 
ing of the system and introduces ZSLC-collisional heat- 
ing and loss processes. But this would only increase (3, 
and we take the observed value as an upper limit for the 
rate constant for FSC and PQNC coUisional loss in this 
state. The scattering length extracted from the simula- 
tion, a « 600 ao, is less reliable because of its sensitivity 
to the trap depth. 

The Pin for Sr ^Pq, which can only reflect PQNC col- 
lisions, is a factor of two larger than /3i„ for optically 
trapped Yb '^P2 [l] and a factor of 5 less than in Sr ^P2. 
The states are sensitive to PQNC and FSC processes. 
ZSLC coUisional loss rates in magnetically trapped 
are one order of magnitude greater |6j], as found in mag- 
netically trapped Ca 

One expects higher PQNC collision rates for •^P2 col- 
lisions than for ^Pq collisions because the process ■^Pj + 
3Pj (bsYSo + [dsGs fSi + AE is allowed for J = 2 
{AE/kB = 1091 /C) and energetically suppressed for 
J = {AE/kB = -581 K). (See Fig. [fl) In light of the 
observed values of /3i„ this implies that PQNC collisions 
make a significant contribution to the inelastic collision 



rate for ultracold ^P2 states. And as predicted @, ZSLC 
collisions should be much more rapid. 

In conclusion, this work has highlighted the important 
role of PQNC collisions in ultracold metastable triplet 
levels in two- valence-electron atoms. It has also provided 
the first measurement of the inelastic and elastic collision 
properties for ultracold atoms in a ^Pq state and an es- 
timate of coUisional properties for •^P2 Sr atoms. These 
levels are of significant current interest in applications 
such as atomic clocks and in fundamental studies in ul- 
tracold atomic physics. The techniques developed here 
demonstrate the production of a cold and dense sample 
of ■^Pq Sr atoms, which is a prerequisite for photoassocia- 
tive spectroscopy of atoms in this state and may allow 
determination of dipole matrix elements and black-body 
radiation shifts that limit the accuracy of optical clocks 
til i22j- 
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